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The large time, small mass, asymptotic behavior of the average mass distribution P{m,t) is 
studied in a d-dimensional system of diffusing aggregating particles for 1 < d < 2. By means of both 
a renormalization group computation as well as a direct re-summation of leading terms in the small 

reaction-rate expansion of the average mass distribution, it is shown that P{m, t) ^ jj I '"^ 1 

for m <^ f^^^, where eKR = e -I- 0{e^) and e = 2 — d. In two dimensions, it is shown that 
P{m,t) ~ iii("iyn(f) £^j. ^ ^ t/ln{t). Numerical simulations in two dimensions supporting the 
analytical results are also presented. 
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I. INTRODUCTION 

Reaction-diffusion systems in low dimensions provide 
an excellent testing ground for developing our under- 
standing of the fluctuation effects in complex systems 
far from equilibrium. A great deal of information, both 
numerical and analytical, has been gathered over the 
past twenty years to show that the evolution of statis- 
tical properties in simple reaction diffusion systems in 
low dimensions is anomalous, in the sense that it does 
not follow the corresponding mean field equations (see 
101 for a short review). In low enough dimensions, in the 
instances where an exact solution is lacking, there are 
no formal methods by which the exponents characteriz- 
ing the different physical quantities may be calculated. 
The renormalization group method (see for a review) 
provides the only systematic way to calculate these ex- 
ponents and thus understand fluctuation-dominated ki- 
netics in reaction-diffusion systems. In this paper, we 
consider the model of irreversible aggregation of diffus- 
ing, massive particles Ai + Aj — > ^i+j, in dimensions 
1 < d < 2, and use the renormalization group method to 
calculate the small mass (m <C f^^^) behavior of P(m, t), 
the mean density of particles of mass m at time t. As we 
explain later in this section, this problem of determining 
P{m, t) requires using the full power of the renormaliza- 
tion group method. This is unlike many other problems 
where just considering the rate equations with a renoma- 
lised reaction rate A is enough to obtain the right answer. 

We now give a more precise deflnition of the cluster- 
cluster aggregation (CCA) model Ai + Aj Ai+j and 
review earlier relevant results. Consider a d-dimensional 
lattice and particles that possess a positive mass. Given 



a configuration of particles on this lattice, the system 
evolves in time via the following microscopic moves: (i) 
With rate D, each particle hops to a nearest neighbor and 
(ii) with rate A, two particles on the same lattice site ag- 
gregate together to form a new particle whose mass is 
the sum of masses of the two constituent particles. As 
time increases, the number of particles decreases due to 
collisions and ultimately when i — s- oo, all particles coag- 
ulate together to form one massive aggregate. However, 
at finite times there is a well defined average mass dis- 
tribution P(m, t) which is of interest to determine. It 
will be shown later (see the text following Eq. (p^)) 
that the large time limit of this model is the same as 
the large A limit. The A — s- oo limit was studied nu- 
merically by Kang and Redner in one and two dimen- 
sions Q. It was shown that P{m,t) has the scaling form 
P{m,t) ~ t~'^f{mt~'^^'^), where d < 2 is the dimension. 
The two exponents are easily determined from the two 
conditions J mP{m)d'm ^ t'^ (mass conservation) and 
/ P{m)dm ~ t~'^^'^ (recurrence of random walks). The 
small mass behavior of P{m, t) can be obtained by know- 
ing the small x behavior of the scaling function f{x). On 
the basis of numerical simulations it was conjectured in 
Q that f{x) ^ j;(2-d)/c; one dimension, the model 
can be solved exactly ||, ||; it was shown that f{x)^x 
or equivalently P(m, t) ~ mt~^^'^ for m <C \/t. The 
one dimensional solution uses the property of ordering 
of particles on a line and is not generalizable to higher 
dimensions. In two dimensions, f{x) was seen numer- 
ically to increase with x for small x Also, in two 
dimensions, the scaling function could be determined in 
the limit of fixed x for t — > oo, where x = mln(i)/t. In 
this case it was shown [|| that P{m,t) = t^^lii^{t)e^^ 
for X < ln^/^(t) and |ln(a;)| < Un(t)|. This result how- 
ever, becomes incorrect in the limit when m is fixed as 
t —^ oo. 
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In this paper, we compute Plm, t) in I < d < 2 in 
the limit t oo, m/rriQ = fixed, where mo is the mass of 
the lightest particle at i = 0. We show that P{m,t) ~ 



2 



and 



e = 2 



ln(r 



for m ^ i'^/^, where ckb = e + O(e^) 
d. In two dimensions, it is shown that 
1^^^ for m <C t/ln(t). These results pro- 



t2 



P(m,t) 

vide a theoretical basis to the results obtained by numer- 
ical methods in Q. 

The CCA model may also be considered to be a special 
case of the more general model in which the aggregation 
kernel is mass dependent. For a review of results on the 
rate equation approach to this problem see ||^, |j . The de- 
pendence of P(m, t) on TO in one dimension in this more 
general model has also been studied In this paper, 
we will restrict ourselves to the aggregation kernel which 
is mass independent; i.e., the rates A and D are indepen- 
dent of mass. When the mass is ignored and A ^ oo, 
the CCA model reduces to the well studied A + A ^ A 
model The CCA model and its variants also find 

application in a large number of physical systems includ- 
ing colloidal suspensions |Tl| , irreversible polymerisation 
[ p^ , aerosols and cloud formation |l2| , river networks 
and coarsening phenomena |l4[ . 

Field theoretic methods have been previously used to 
study complex systems far from equilibrium (see for 
a review). We briefly review results relevant to reaction- 
diffusion systems. In some earlier works the ef- 
fective reaction rate and the decay exponent of the aver- 
age particle density were computed for the A + A^ ^(0) 
model. The renormalization group study of the same 
model with sources was done in In the sys- 
tematic renormalization group procedure for the compu- 
tation of average density and density-density correlation 
function in kA % reaction was developed. In |Q renor- 
malization group analysis of A + B ^ % reaction in d > 2 
was used to study the effects of initial fluctuations on 
the late time decay of particle densities. The renormal- 
ization group technology developed by Peliti, Lee and 
Cardy [0 ^ was used to compute the average mass 
distribution of clusters in the CCA model in the inter- 
mediate mass range in 

It turns out, however, that as far as the study of scal- 
ing properties of one-point correlation functions in most 
reaction-diffusion systems is concerned, renormalization 
group is not a vital tool. Consider, for example, a sin- 
gle species annihilation model A + A 0. Once the 
renormalization of effective reaction rate is understood, 
the correct density decay exponent can be obtained from 
simple dimensional arguments [ p^ . Alternatively, one 
can use simple random walk arguments or use the Smolu- 
chowski approximation (we refer to the case in which 
the reaction rate is replaced by a time dependent reac- 
tion rate as the Smoluchowski approximation) to 
obtain the correct values of decay exponents. A renor- 
malized mean field theory or, alternatively, a version of 
Smoluchowski theory can also be used to compute the 
average mass distribution in cluster-cluster aggregation 
for intermediate masses [ p2[ . 

In the CCA model considered in this paper, it turns 
out that the stochastic field P{m,x,t), describing the 



continuous limit of the local mass distribution, has a 
non-zero anomalous dimension in d < 2. The scal- 
ing exponent governing the dependence of the average 
mass distribution P{m, t) on mass is proportional to the 
anomalous dimension of the operator corresponding to 
the local mass distribution. As explained in Section || 
any approximation scheme which disregards this anomaly 
(such as the Smoluchowski approximation) predicts that 
Pijn, t) ~ when m <^ f^l"^ , for any dimension. This is 
in contradiction with both numerical results as well as 
the exact result in one dimension ^|. The full power 
of renormalization group analysis has thus to be brought 
to bear in order to compute the anomalous dimension of 
P(to, X, t) in the form of the e = 2 — d expansion. Since 
the peculiarity in the small mass distribution of cluster- 
cluster aggregation is first discussed in we refer to this 
as the Kang-Redner anomaly. 

In addition to the calculation of P(m,t), we found 
many similarities between the problem of cluster-cluster 
aggregation and the problem of weak turbulence (see |2|] 
for a review on weak turbulence). We elaborate on this 
connection in Sec. VI, where the Kang-Redner anomaly 
is interpreted as a violation of the Kolmogorov (constant 
flux) spectrum of particles in mass space due to strong 
flux fluctuations. 

The rest of the paper is organized as follows. Section || 
contains a discussion of the stochastic integro differential 
equation satisfied by P{m,x,t). The mean field results, 
as well as the reasons for their failure in low dimensions 



are also included. In Sec. Ill, we analyze the large time 
asymptotic behavior of P(to, t) in d < 2 using the renor- 
malization group method. We do the same for d = 2 in 
Sec.0. In Sec.0, we redcrive the results of Sec. 03 and 
IV using explicit re-summation of all diagrams giving the 
leading contribution to the average mass density in the 
limit of large time. Reasons for the failure of the Smolu- 
chowski theory then become more transparent. In Sec. VI 
we elaborate on the connections with weak turbulence. 
Finally, we conclude with a summary and discussion of 
open problems in Sec. VII. 



II. CONSTANT KERNEL CLUSTER-CLUSTER 
AGGREGATION AND MEAN FIELD ANALYSIS. 

The problem of computing density correlation func- 
tions in d-dimensional stochastic processes can be refor- 
mulated as an effective equilibrium problem in (d + 1)- 
dimensions with the help of Doi-Zeldovich-Ovchinnikov 
trick One can then attempt to solve the problem 

using the powerful methods of statistical field theory, in 
particular those based on renormalization group ideas. 
Starting from the lattice version of the CCA model, 
we would like to derive the corresponding field theory 
and the Langevin equation obeyed by P{m, x,t), where 
P{m, X, t)dmdV is the number of particles with masses in 
the interval [m,m-|-dTO] in the volume dV. It was shown 
in j6| , that the problem of finding all correlation functions 



3 



of the local mass distribution is equivalent to finding all 
moments of the following stochastic integro-differential 
equation (stochastic Smoluchowski equation): 



d_ 

dt 



DAj P{m,x,t) ^XP*P 
- 2XN{x, t)P{m, X, t) + 2^^/X^{x, t)P{m, x, t), 



(1) 

where P * P = dm' P{m — m' , x,t)P{m' ,x,t) is a 
convolution term, ^{x, t) is Gaussian white noise and 
N{x,t) = d'mP{m,x,t), is the local particle den- 
sity. We are interested in P{m,t) = {P{m,x,t)), where 
(. . .) denotes averaging with respect to the noise ^. If 
the initial number of particles of different masses at 
different lattice sites are independent Poisson random 
variables parameterized by initial average mass distri- 
bution PQ{m), then the initial condition to be supplied 
with Eq. (|l|) is P{x,m,0) ~ PQ{rn). It is easy to 
check that Eq. (Q) conserves the average mass density 
P = J mP{m)dm. 

As particles aggregate the typical mass grows in time 
as t''/^. If we are interested in small masses, we need 
to consider masses smaller than the typical mass, m ^ 
p{Dty/^. This can be achieved by considering m/rriQ 
to be fixed as i oo, where mo is the smallest mass 
at t = 0. In this case, the first term on the right hand 
side of Eq. (Q) is almost surely small compared to the 
other terms. Consequently, the small mass behavior of 
the local mass distribution is described by the following 
system of non-linear stochastic partial differential equa- 
tions (SPDE): 

^-m^ P{m,x,t) 



= -2\N{x,t)P{m,x,t) 
+ 2iVX^{x,t)P{m,x,t), (2) 
= -XN^{x,t) 

+ 2iVX^{x,t)N{x,t). (3) 
Equations and demonstrate an interesting con- 



§-^--DA\Nix,t) 



ncction between this model and A + A 



A model. 

Stochastic field P can be identified with S^, where A^o 



is the initial density. Differentiating Eq. (|3|) with respect 
to No, and setting P = we obtain Eq. (^). 

We are interested in the behavior of P(m, t) in the 
limit of fixed m and t —^ oo. We can then identify the 
particles with this fixed mass as B kind of particles and 
the remaining particles as A kind of particles. Then, 
clearly, the study of Eqs. (^) and (^) is equivalent to the 
study of a two species reaction 



A + A 
A + B 



A, 

Inert, 



(4) 



in the limit when concentration of B particles is much 
smaller than that of A particles. This two species prob- 
lem has been studied in d = 1 for arbitrary diffusion rates 



p5[ . Specializing results of this paper to our case, we find 
that P{m, t) ~ t~^/'^ for t — > oo. Assuming that the large 
time asymptotics of P{m, t) is universal, we can restore 
the m-dependence using dimensional analysis, to obtain 



P{m,t) = C 



m 



(5) 



where C is a constant, and p is the average mass density. 
Equation ^ matches with the exact results obtained for 
the CCA model in one dimension ||. 

However, no exact solutions are available for dimen- 
sions d > 1. In the rest of the paper we will be analyzing 
Eqs. (|^) and (^) in d > 1 using the dynamical renor- 
malization group method. We will show that for small 
masses and 1 < d < 2, 



P{m,t) 



1 ( 



piDty \p{DtY/'^ 
where ckb. = e + 0{e^) and e = 2 — rf. li d^2 



m^p{Dtf/\ (6) 



P{m,t) 



■ln(t/io)ln 



for pDta ^ m <t: pDtlii(t/to 



1 



\n{t/to 



(7) 



Here, ~ A'^/D, where A is the lattice spacing. 

Before doing the renormalization group analysis, let us 
analyze Eqs. and (^) in the mean field (weak coupling) 
limit. Neglecting stochastic terms in the right hand side 
of Eqs. (g) and (^ and solving the resulting system of 
ordinary differential equations, we obtain 



NMF{t) = 
pMFit) = 



Nn 



1 + NoXt ' 
(1 + NoXty 



(8) 
(9) 



Thus, at large times, P{t) ^ t , given that mean 
field theory is applicable. Relative corrections to the 



At 



There- 



mean field result are of the order (t) — ■ 
fore mean field theory is asymptotically exact in d > 2 
p6t and diverges with t if d < 2. Re-summation of the 
most divergent terms in the weak coupling expansion of 
P around the classical solution is required and can be 
performed in the case at hand using the formalism of 
renormalization group. The details of the computation 
are given in the next Sec. III. Here, we would like to 
demonstrate that accounting for renormalization of the 
coupling constant alone doesn't yield the correct decay 
exponent as mentioned earlier. To the leading order in 
e = 2 — d, renormalization of the effective reaction rate 
reduces to replacing A in Eqs. (|^) and (^) with renormal- 
ized value = f{e)t~^/'^ and omitting stochastic terms 
(renormalized mean field approximation). Eliminating 
N{t) from the resulting system of ordinary differential 
equations one finds the following equation for P{m, t): 



dP 
'dt 



P 



(10) 



N 



N 




(A) 



■X = 



-X + 



•ft 




(B) X = 



-X + 



FIG. 1: Propagators and vertices of the effective theory, 

Eq. 0f (C) 



This implies, that P{m, t) ^ mPt^'^. In other words, 
P(m, t) scales with time according to its physical dimen- 
sion. As a result it does not depend on mass. As we will 
show in Sec. Ill, arguments leading to this conclusion are 
incorrect, as they disregard the possibility of anomalous 
dimension of the stochastic field P. 



(D) 



+ 



FIG. 2: Diagrammatic form of mean field equations for: (A) 
Average particle density; (B) Average mass distribution; (C) 
A'^A''-response function; (D) PP-response function. 



III. RENORMALIZATION GROUP ANALYSIS 
OF STOCHASTIC AGGREGATION. 

The average mass distribution P{t) and average parti- 
cle density N(t) admit functional integral representations 
which can be obtained by applying the Martin-Siggia- 
Rose procedure J 27j to Eqs. (||) and (||) (equivalently see 
Eqs. (2)-(4) of 1 1). We then obtain 



VN{x', t')VN{x', t')VP{x', t')VP{x', t') 



X 0{t) 



-S^ff[N,N,P,P] 



(11) 



where 0{t) — N{x, t) or P{x, t) and 

5gff= I d'^xdT[N{N - DAN) + P{P - DAP) 
Jo 

+ X{NN^+ 2PPN + N^N^+ 2NPNP + P^P^)] 

(12) 

is the effective action functional. Perturbative expan- 
sions of N{t), P{m, t) in powers of A can now be obtained 
in the standard way, see for example ||2^ . Feynman rules 
for constructing terms of these expansions are summa- 
rized in Fig. |l|. Due to the non-renormalization of the 
diffusion rate as well as the average mass density in the 
field theory Eq. ( |l^), in all that follows, we use units in 
which p — D ~ 1. 

The average mass distribution P(m, t) is formally given 
by the sum of all diagrams built out of blocks presented 
in Fig. |l| with one outgoing punctuated line (P-line) . The 
contribution from each diagram is a function of AiVo^ and 
'bare' dimensionless reaction rate go{t) = ^jf^- Unless 
we are interested in the small time expansion of P{m, t), 
XNot cannot be treated as a small parameter. Therefore, 



the contributions of all diagrams proportional to a given 
power of go{t) and various powers of AA^o^ have to be 
summed up. A simple combinatorial argument (see 
for details) shows that the contribution of a diagram with 
n loops is proportional to gQ^t)". In the weak coupling 
regime the main contribution to P(rn, t) and N(t) is, 
therefore, given by the sum of all tree diagrams, the first 
correction comes from summing all one-loop diagrams 
and so on. It turns out ^ that the sum of all tree 
diagrams gives the mean field answers Eqs. (H) and (||). 
The expansion in powers of go{t) is therefore the standard 
loop expansion around solutions of mean field equations. 
It is obvious that in d < 2 the loop expansion is not 
very useful at large times as limi^oo5o(i) = oo. Fortu- 
nately, the value of dimensionless reaction rate properly 
renormalized to account for the build-up of inter particle 
correlations turns out to be of order e = 2 — d for large 
times. This allows one to convert the loop expansion 
into an e-expansion using perturbative renormalization 
group method. Such an expansion works well for large 
times and will therefore yield all the information we need 
about the behavior of the average mass distribution in 
the strongly fluctuating regime. 

In computing loop corrections to any order, there are 
generally an infinite number of diagrams to sum. How- 
ever these diagrams can be resummed in part if one in- 
troduces classical response functions and G^^. Re- 
sponse function G^^ [G^f) is equal to the sum of all 
tree diagrams with one outgoing and one ingoing line 
of types N (P). As was already mentioned, mean field 
densities Eqs. ^ and (^) are also equal to the sums of 
tree diagrams with one outgoing N or P line correspond- 
ingly. One then simply has to associate incoming lines 
with mean field densities and propagator lines with mean 
field response functions. Integral equations satisfied by 
classical densities and response functions are presented 
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FIG. 3: One-loop corrections to the mean field answer for 
average mass distribution. 



in diagrammatic form on Fig. |2[ 

Solutions of Eqs. Fig. ||(A) and Fig. ||(B) coincide with 
Eqs. (|) and (||), as they should. Equations Fig. ||(C) and 
Fig. |2(D) can also be solved with the result 



G 



G^l {x2,t2]Xi,ti) = 
I 

Gq{x2- Xi,t2-ti), (13) 



Noit2 

where Go is the Green's function of the linear diffusion 
equation. 

Using the notion of mean field response functions and 
densities one can easily classify all one-loop diagrams 
contributing to average mass distribution. The result 
is presented in Fig. ^. A quick check shows that analyt- 
ical expressions corresponding to diagrams Fig. ^(i) and 
Fig. |3|(iii) containing primitive loops diverge in d > 2, 
which is consistent with that fact that upper critical di- 
mension of the effective theory Eq. ( |T^ ) is two. Comput- 
ing the relevant integrals in d ~ 2 — e dimensions we find 
the following one-loop expression for the average mass 
distribution: 




1 



;F{e) 



where F(e) 



8^)d/2 Af(2-E/2)- 

2-loop corrections, (14) 



l-e/2 

Equation ( |14|) can be used to determine the large time 
asymptotics of P{t) in the following way. The exact aver- 
age mass distribution satisfies the Callan-Symanzik equa- 
tion which we will derive below. Coefficients of this equa- 
tion depend on the law of renormalization of all relevant 
couplings of the theory Eq. (12). One relevant coupling is 
the effective reaction rate. Its renormalization is known 



|r7[ . Below we will show that the only other relevant 
coupling is the initial mass distribution Pq. We will de- 
termine it's renormalization law with one-loop precision 



demanding that, the expression Eq. (|lj) be non-singular 
in the limit e ^ when expressed in terms of renor- 
malized coupling constant and renormalized initial mass 
distribution. This will determine coefficients of Callan- 
Symanzik equation up to terms of order gj^, where 
is the renormalized dimensionless reaction rate. Solving 
this equation we will be able to compute the decay expo- 
nent of P(t) up to terms of order e. 

Dimensional analysis of effective vertex parts of the 
theory Eq. ( [l^ ) shows that the only relevant bulk cou- 
pling is the effective reaction rate. Its relevance is due to 
the recurrence of random walks in d < 2. Reaction rate 
renormalization accounts for all fluctuation effects in the 
A + A model. 

However, the CCA model is more complicated and in- 
teresting due to the presence of 6oun(iar?/-relevant cou- 
plings. To identify them, we use the following version 
of dimensional analysis. Boundary coupling constants 
correspond to vertices with no incoming lines. As PP- 
interactions can be neglected in the problem at hand, 
we are interested in boundary vertices with at most one 
P-line. Assume for simplicity that d = 2 (critical dimen- 
sion). Assume also that the initial density A^o = oo. This 
assumption is justified if one is interested in the large 
time behavior of correlation functions in aggregation, as 
Nq flows to infinity under renormalization group transfor- 
mation to increasingly larger time scales, see jist for more 
details. Let Ta^pit), where a = 0, 1; = 0, 1, 2, . . . be 
the simultaneous Green's function of the theory Eq. ( p^ ) 
with a P-lines and f3 N-lines with all external momenta 
equal to 0. Using Eq. (p^), one can express Ga,[3{t) in 
terms of corresponding vertex part Va^pit) as follows: 



G 



a+l3 



drr2("+«F„_^(r). (15) 



As physical dimension of Ga.f}{t) is (Length)~^^~'^'^ , 
Va,i3{t) ^-/3-2q-i_ result, Eq. ( p^ ) converges at 

small times for any a if /? > 0. If /3 = and a = 1, 
Eq. ([l^) diverges logarithmically. This divergence can be 
regularised using a small-time cut-off and leads to 
the renormalization of the initial mass distribution Pq. 
The latter plays a role of the coupling associated with 

1^1,0- 

As a result of the discussed divergence, diagrams in- 
volving Vi^o grow with time faster compared with dia- 
grams with the same number of loops but with no sub- 
diagrams contributing to renormalization of Pq. Conse- 
quently, diagrams belonging to the former class have to 
be re-summed exactly in order to obtain the correct large 
time behavior of the average mass distribution. 

We conclude that fluctuations in stochastic aggrega- 
tion lead to two effects: renormalization of the effec- 
tive reaction rate and renormalization of the initial mass 
distribution. It follows that renormalization of these 
two couplings regularizes the perturbative expansion of 
P{m, t) to all orders. 

As it turns out, the renormalization of Pq is solely re- 
sponsible for Kang-Redner anomaly. In Sec. ^ we will 
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analyze initial density renormalization by explicitly re- 
summing small time divergences in the perturbative ex- 
pansion for P{in,t). Now, we will show how it appears 
formally within the framework of perturbative renormal- 
ization group method. We follow dynamical renormal- 
ization group procedure described in ||l5|] , part I. For a 
review of results concerning coupling constant renormal- 
ization, see ibid, part II. 

Let us fix a reference time io > 0. Expression ( p^ ) 
evaluated at is to be made finite by absorbing the 
divergences appearing as e ^ into the renormalization 
of the reaction rate and the initial density. 

Let go = -^^0^^ be the dimensionless 'bare' reaction 
rate. As has been shown in [ pT| , renormalized reaction 
rate gn is related to go via the following exact formula: 



gR 



90 



1 



90 
9* 



(16) 



Here g* — 27re(l + 0(e)) is the non-trivial fixed point 
of the renormalization group flow in the space of effec- 
tive couplings of Eq. (|lj). Recall that go ^ A. Hence 
Ihnx^oo gR = g* ■ It follows from the Callan-Symanzik 
equation that large time behavior of the P(to, t) is also 
determined by the fixed point value of the effective re- 
action rate g*. We therefore conclude that the limits 
i > 0, A — > oo and A > 0, t — > oo of Kang-Redner model 
belong to the same universality class as claimed in the 
introduction. 

Solving Eq. ( p^ ) with respect to go, expanding the re- 
sult in powers of gji and substituting the expansion into 
Eq. (^) we obtain the average mass distribution at time 
to as a power series in gji: 



Pito) = 



Po 



NS9%4 



1 



1 

7 



'l + Oig*))gR + 0{gl) 



.(17) 



As expected, the order-g^j term in Eq. ( |17| ) is still singular 
at e = 0. To cancel the remaining divergences we have 
to introduce renormalized initial mass distribution Pr: 



Z{gR,to,e)Po, 



(18) 



where Z is a power series in gji with coefhcients chosen 
in such a way, that the average mass density 

P{t,gR,PR,to) = Z{gR,to,e)P{t,Xo,Po,e). (19) 

when expressed in terms of PR,gii, is non-singular at 
e = 0. Substituting Eq. ^ into Eq. ^ we find that 



Z=l+'J^ + 0{gj,), 



(20) 



in order for P{t) to be non-singular at 1-loop level. 

Now we can derive the Callan-Symanzik equation. The 
fact that P{t, Ao, Po, e) does not depend on the reference 
time to leads to the following equation for P{Pr): 

to^{Z-^P{Pn)) = 0. 

Oto 



Noticing that P{t) = to'^jfi^(t/to, gn), where $ is a 
dimensionless function, one can convert the above condi- 
tion into Callan-Symanzik equation for P(t): 

+ ll^(aR)^+d-^jigB)^P{t,gR,PR,to) = 0, 

(21) 

where P{gR) — —2-^gR is the beta function of the theory 



Eq. (y) and 

i{9r) - 



1 

Z' 



d_ 

'dto 



-2-to^Z = - 



9R_ 
27r 



+ 0{gl,gRe), (22) 



is the gamma function. 

It is well known, that at large times and in d < dc — 2, 
solutions to Eq. ( pi[ ) are governed by non-trivial fixed 
points (zeroes) of beta-function. Differentiating Eq. ( [l^ ) 
with respect to to, one finds that P{gR,) — gR{gR — g*)- 
Hence, (3{gR) has a unique non-trivial fixed point gR — 
g*. It follows from Callan-Symanzik equation (Eq. (|2l|)) 
that P(t) ~ t"''*, t oo, where 



d* = d- ^j{g* 



(23) 



We see that scaling dimension of P{m, t) differs from it's 
physical dimension by a term proportional to the value 
of 7-function at the fixed point. This term is called the 
anomalous dimension of the field P. The physical reason 
for the presence of anomalous dimension is the renor- 
malization of the initial mass distribution by small time 
fluctuations. 

As g* e, the substitution of gR expansion of 7 into 
Eq. ( P3|) yields the e-expansion for d* . Using Eq. (^) we 
flnd that 



d* =d+^e + 0{e^). 



(24) 



Finally, one can restore the m-dependence of P{fn, t) us- 
ing dimensional argument. The result is 



1 fm^Y"'' 

P^„,.t)^A^.)-,(-^) . 



(25) 



where ckr — e + O(e^) is equal to twice the anomalous 
dimension of the stochastic field P. 

Note that in d = 1, ckr = 1 and P(m, t) ^ m} , which 
coincides with the exact answer, However, at the 
moment we do not have reasons to believe, that higher 
order e-corrections to our answer for Ckr vanish iden- 
tically for any < e < 1. Equation ( p4|) implies that 
anomalous dimension of P vanishes aX dc — 2. Yet, it fol- 
lows from general theory of perturbative renormalization 
group that traces of anomalous dependence of P(rn, t) 
on mass in d < 2 must be present at the critical dimen- 
sion as well. We analyze Kang-Redner anomaly in d = 2 
in Sec. IV and compare renormalization group predic- 
tions with conclusions of direct numerical simulations of 
the system of diffusing-aggregating point particles on the 
two-dimensional lattice. 
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FIG. 4: The variation of P{m, t)t^ with t is shown for m = 1, 2 
and 4 on a semi-log scale. The variation is linear implying that 
P{m,t) — c(m) ln(t)/f^, where c(m) is some mass dependent 
function. 



FIG. 5; The small m behavior of P(m,t) / P(l,t) is shown on 
a semi-log plot. The variation is linear implying that c(m) 
ln(m) with c(m) as defined in caption of Fig. 



IV. KANG-REDNER ANOMALY IN TWO 
DIMENSIONS. 

It is well known that anomalies 'mdc—e dimensions lead 
to logarithmic corrections to mean field theory answers in 
d — dc- Kang-Redner anomaly is not an exception. Small 
mass behavior of the average mass distribution in two 
dimensions can be easily obtained by solving the Callan- 
Symanzik equation (see Eq. (|2l|)). Note that 

Pig) \d=2 = -^9^, 

7(5) U=2 = -^<?(l + 0(.g2)). (26) 

Solving Eq. (^l|) with coefficients given by Eq. (^6|) and 
the initial condition 

Pito) = Const^(l + O(gR)), (27) 
Sr'-o 

produced by the mean field theory, we find that 

P{t) = Const ^"^^f°V + O(l/ln(t/to)). (28) 
t^ 

Such a behavior of P{m^t) \m=fixed in two dimensions 
was originally seen by Kang and Redner in numerical 
simulations . We have now shown that Eq. (^8|) can be 
obtained as a result of systematic renormalization group 
computation. 

Analyzing their model in two dimensions, Kang and 
Redner noticed that P(m, t) is a slowly varying (increas- 
ing, as suggested by Fig. || of [||) function of mass. We 
can now quantify this observation by computing the mass 
dependence of the average mass distribution. In d = dc 
direct dimensional arguments cannot be used to restore 
the mass dependence of Eq. (^) due to the explicit de- 
pendence of P(to, t) on the lattice cut-off. Instead let us 
analyze Eq. (p5|) for P(m, t) in the limit e ^ -1-0. Near 



e = 0, the amplitude A(e) ~ e~^. Expanding the right 
hand side of Eq. ( psf ) in powers of e and using the fact 
that P{m, t; e) is non-singular at e = as a function of 
renormalized parameters, we find that 



P(m, t) = - [Ci ( ln(t/io))' + C2 ln(t/to) In(TOAo) 

+C3{\n{m/h))A (1 + O{l/\n{t/to))).{29) 



We are interested in P{m,t) in the limit t — s- 00, m = 
fixed. Time dependence of P(m,t) is given by Eq. ([2^). 
Therefore, coefficient Ci in Eq. ( p^ ) is 0. Hence 

P(m, t) = Const ^"(^/^°\t'^("^/^°) (l + 0(1/ ln(t/to))) . 

^ (30) 
Note that Eq. (gg) is valid for m < M{t), where 
M{t) ^ t/ ln(f/fo) is the typical mass. If to ~ M{t), then 
P{m,t) ^ {ln{t / to) / 1)'^ , which coincides with the answer 
for P(m,t) Im^M(t) ij^ the intermediate mass range 

To check our analytical results Eq. (|2^) and Eq. ( |30| ) 
we did a numerical simulation of the system on a two- 
dimensional lattice. For the sake of computational effe- 
ciency, we chose to work in the limit A — *■ cxd so that a 
lattice site can hold utmost one particle. The lattice was 
chosen to be a square lattice of size 3000 x 3000 with 
periodic boundary conditions. In Fig. ^ we show the 
variation of P{m, t) with t for fixed small to. It is seen 
that P{m,t) — c(to) ln(i)/i^, where c(to) is some func- 
tion of the mass to. This is in excellent agreement with 
Eq. (p8|). To determine c(to), we studied the variation 
of P{m,t)/P{l,t) with the mass to (see Fig. |]). We see 
that c(to) ~ In(TO), thus confirming Eq. (pO|). 
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V. KANG-REDNER ANOMALY VIA EXPLICIT 
RESUMMATION OF BOUNDARY 
DIVERGENCES 

In this section we will derive Eq. (|2^) without using 
the formalism of renormalization group. Instead, we will 
identify the principal set of diagrams contributing to the 
the large time limit of P(m, t) and derive a simple integral 
equation satisfied by the sum of these diagrams. 

Additional divergences in the terms of the perturbative 
expansion of the average mass distribution discussed in 
the previous Sec. IV are actually due to the quadratic 
singularity of the mean field expression Eq. (|9|) for the 
average mass distribution at A'o = 00. t = 0. Let us 
illustrate this statement on the example of diagram (i) of 
Fig. 1^. Assuming for simplicity that d ~ 2 and denoting 
ultraviolet cut-off by Ai, we find that the most divergent 
contribution coming from this diagram is 



27rA2Ar2i2 



In 



t 

At 



t 

dr- 



XNo 



1 + AA^or 



[1+0(1/ XNot)). 

(31) 

Setting A^o = 00 under the sign of integration is clearly 
impossible, as the resulting integral would diverge. Di- 
rect computation of Eq. (pll) shows that 



In 



^ ) HXNotKl 



0{l/XNot)). 



We see that lu-j is proportional to \n{t)^, not 



(32) 
ln(t) as 

one would have deduced from the simple loop count- 
ing. The additional singularity is due to the fact that 
Pmf |Afo=oo~ i"^- (Recall that Nmf ko=oo~ t^^- As a 
result divergences at t = are absent in the field theory 
of A + A ^ reaction.) The additional boundary singu- 
lar terms are generally present in all Feynman integrals 
contributing to P{m, t) and have to be re-summed to all 
orders of perturbation theory to get the correct large time 
asymptotics of the average mass distribution. 

Let n(i2, ^i) be the exact zero momentum vertex func- 
tion of type PP - the sum of all one-particle irreducible 
diagrams contributing to PP response function divided 
by the propagators corresponding to external lines. The 
Schwinger-Dyson equation for the exact average mass 
distribution reads: 



P{t)=PMF{t) 



1 
^2 



dU 



dhn{t2,h)Pih). (33) 



Equation ( p3| ) is most easily derived using the formalism 
of Feynman diagrams. Its diagrammatic representation 
is given in Fig. ^ (I). 

We know that in two dimensions P{t) ~ t~^. There- 
fore, near d = 2, the function ri{t) = P(t)Xt^ is slowly 
varying. The following simplified version of Eq. (33) is 
therefore valid to the leading order in e: 



Tj{t) = 1+ / dmitiMh), 



(34) 



X -X 

-I- -^^ — B — ■ p azDoooA^ 

t2 \ 7l1 A\ 



12\ / t1 



12\ / t1 



FIG. 6: (I) Schwinger-Dyson Equation for n(t). (II) Pertur- 
bative expansion of the polarization operator: (IIA)- in the 
number of loops; (IIB)- in the order of exact vertex parts. 



where 



n(<) = t2 r^n(i,ti) 

JO ^\ 



(35) 



Differentiating Eq. ( p4D with respect to time, we find that 
r\(t') satisfies the following differential equation: 



dri 
'dt 



(t) = mvit)- 



(36) 



It turns out (see below), that at large times 



0(e2))(l + 0(l/t^/2)). (37) 



Solving Eq. (||) with U{t) given by Eq. 1^), we find 
that P{t) = ri{t)/e - i(-2+£/2+o(e=^))^ rpj^^g^ Yia.Ye 
been able to rederive Eq. ( p^ ) without any reference to 
renormalization group. 

Before we turn to the derivation of Eq. (|37|), we 
would like to demonstrate, that solving integral equation 
Eq. ( p4[) is indeed equivalent to summing leading bound- 
ary singular terms in the perturbative expansion of P{t) 
to all orders. All divergences at i = 0, iVo = 00 can be 
regularised by setting the lower limit of integration in the 
integral in the right hand side of Eq. (^^ equal to . 
Resulting equation, 



77(<) 



dti 



(38) 



can be solved using the method of consecutive approxi- 
mations. Setting ri{t) = ?7o(i) + Vi(t) + %(i) + • • ■ and 
treating the integral in the right hand side of Eq. ( p^ ) 
perturbatively, we find that r]o{t) ~ 1 and 



Vn{t) ^ ^_{^HXNot)y\ n 



1,2, 



(39) 
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Note that rji (t) is equal to 1-loop boundary singular term, 
772(0 " 2-loop boundary singular term and so on. Sum- 
ming the resulting series for j] one confirms our main 
result that r]{t) ~ f^^"^. One loop-boundary singularity 
leads to mass-dependent logarithmic corrections to the 
average mass distribution. This is most easily seen by 
expanding Eq. ( |25| ) in powers of e. These corrections, 
which were originally observed in |^ , motivated to a cer- 
tain extent the present investigation. 

The knowledge of the first term in the e-expansion of 
decay exponent of rjit) leads to a good approximation of 
the average mass distribution , given that eln(AA'ot) ~ 
1, but In(AA'ot) <C 1. If the latter condition breaks 
down, we need to know the expansion of the polarization 
operator Il{t) to second order in e and modify Eq. (jsj) by 
including corrections proportional to the derivative r]{t) 
(which account for the fact that rj{t) is not constant). 

We will now derive Eq. (|3^) by computing Tl{t) to the 
first order in e using 1-loop diagrams and verifying that 
higher- loop diagrams lead to terms of order and higher. 
If mean field theory for the average mass distribution 
was exact, operator 11 would have been identically equal 
to zero. Consequently only loop diagrams contribute to 
Il(t2,ti). The one-loop diagrams are shown in Fig. ^ 
(IIA) . The computation of corresponding Feynman inte- 
grals is straightforward. After integration with respect 
to ti we find their respective contributions to Ii(t): 



Adding together Eqs. (|o|), (Bll) and (||), we find that 



n(.o(0 = 



-8Xt'/^ (1 -I- 2e) 
(87r)'^/2t e(H- e/2)2 ' 

(87r)'i/2te(l + e/2)' 

(87r)'i/2t (l + e/2)(2 + e/2)' 



(40) 
(41) 
(42) 



Note that individual contribution from diagrams (i) 
and (ii) are 1/e-times bigger than the contribution from 
diagram (iii), which does not contain primitive loops. 
Yet, terms of order e^^ cancel upon adding Il(^ij{t) and 
(t) leaving terms of order up to - same as the lead- 
ing order of terms of Tl(^iii) (t). This cancellation explains 
why we had to account for an apparently subleading con- 
tribution of diagram (iii) to the perturbative expansion 
of n(i). Such a cancellation is not accidental and hap- 
pens at all orders of loop expansion: diagrams (i) and 
(ii) can be interpreted as first two terms in e-expansion 
of the first term in the cumulant expansion of 11(^2,^1), 
see Fig. ^ (HB). This term corresponds to an exact total 
particle density N{t) connected to the PP-propagator 
via an exact PTVP-vertex. The vertex is of the order of 
the fixed point coupling g* ^ e, while exact density is 
of the order e~^, |1^. Therefore the contribution of the 
term in question to 11(^2,^1) is of order 1, which is re- 
flected in cancellation of the terms of the lowest order in 
e in every term of its loop expansion. 



m 



2At^/2 



(1 + 0{e)) + 2-loop corrections. (43) 



Large-t behavior of Tl{t) can be obtained from Eq. (|42 
by replacing the 'bare' reaction rate A with renormalized 
reaction rate Ai?(f) ^ 2net~'^^^,t — > 00. The result does 
indeed coincide with Eq. (|3^). 

It remains to verify that two- and higher loop diagrams 
contribute only to higher order terms in the e-expansion 
of n(i). Order-e~2 contributions from diagrams con- 
taining only primitive loops cancel as explained above. 
Order-e~^ contributions from diagrams with two primi- 
tive loops and two-loop diagrams with one primitive loop 
are accounted for by one-loop renormalization of cou- 
pling constant in one- loop diagrams. Hence non-trivial 
corrections to polarization operator come only from two- 
loops diagrams containing no primitive loops and non- 
singular parts of all other two-loop diagrams. Simple 
counting shows that contribution from such diagrams to 

Tl{t) is proportional to '•''^"ff_*2)f ~ t- Similar argu- 
ment shows that n-loop diagrams contribute to Tl{t) at 
the order e" only. 

Now it is very easy to characterize the class of dia- 
grams giving the leading contribution to the e-expansion 
of n(i). The statistics of N{t) is strongly non-Gaussian. 
Yet, the main contribution to the polarization operator 
n(i) comes from the diagrams proportional to the first 
and the second cumulants of the stochastic field N{x, t) 
only. Non-Gaussian effects are due to the fact that these 
cumulants are connected to PP-propagator via exact ver- 
tices. 

It is also possible to derive formula Eq. (|2^) with- 
out using the formalism of renormalization group. In- 
stead, one can solve equation Eq. ( |36| ) directly using 
one-loop expression for polarization operator Hit) in two 
dimensions. The latter can be obtained by setting in 
Eq. ( ^3| ) d = 2 and replacing the bare reaction rate 
A with renormalized reaction rate in two dimensions: 
Mt) = 15(77*^(1 + O(l/ln(t/to))). The presented ex- 
pression for Afl(t) is easy to compute by explicit re- 
summation of all diagrams contributing to renormaliza- 
tion of the bare reaction rate, see ||l^ for details. The 
resulting equation for ri{t) = t^P{t) is 



1 



tHt/to) 



(44) 



The solution is r]{t) ~ ln(i/io). Correspondingly, P{t) 



, which coincides with Eq. (p8D 



VI. KANG-REDNER ANOMALY AND 
CORRECTIONS TO KOLMOGOROV PARTICLE 
SPECTRUM. 

Now we will show that the mean field result P(m, t) ^ 



m 



can be interpreted as a constant flux (Kolmogorov) 
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solution of the Smoluchowski equation. We will then in- 
terpret Kang-Redner anomaly as a breakdown of Kol- 
mogorov scaling due to strong flux fluctuations develop- 
ing at large times. 

Averaging Eqs. (||) and (^) with respect to noise one 
gets the following relation between one- and two point 
mass distribution functions: 



d{P) 
dt 

d{N) 
dt 



= -X{PN), 



(45) 
(46) 



Let us look for solutions of Eq. (|46|) having the form 



P{m,t) 



m 

M{t) 



(47) 



where ^ = Recall that M{t) = ^ is the typ- 

ical mass. The new dependent variable J has a simple 
physical meaning: dfi'J{^', t) is the average number of 
particles with masses less than M(t)fi contained in the 
volume N^^{t). 

At times much less than tc = (A)"^/*^ relative fluctua- 
tions of local density are small. As a result, mean field 
theory is applicable and (JN) « {J){N). As a result, 
equation Eq. (ksf) simplifies to 



-(-) 



dJ_ 



(48) 



Therefore, ^ is a locally conserved quantity with flux 
equal to (—J). Note that J > 0. Therefore the cascade 
of -jj is inverse in the terminology of turbulence: its flux 
is directed towards the small masses. We see that self- 
similar solutions of Eq. ( ^ ) correspond to constant flux 
solutions of Eq. (^8|). The latter is just J = Const. Con- 
stant flux solutions of kinetic equations are called Kol- 
mogorov solutions in the theory of weak turbulence, see 
p3[ for details. We therefore conclude that in the mean 
field approximation 



(P)(m,<) 



Njt) 
M{t) 



(49) 



where ekoim = is the exponent, which determines the 
Kolmogorov scaling of the average mass distribution. 
Note, that the fiux J also has a meaning of dimensionless 
particle density (the number of particles in the volume 
, which is an obvious integral of motion. 
The fact that J = const, means that particles are 
equipartitioned between system's degrees of freedom and 
particles' flux is identically equal to zero. The character- 
istic feature of the state Eq. ( p9| ) of our system is there- 
fore the presence of non-zero constant flux of one integral 
of motion and equipartition of the other. Similar kind of 
behavior has been observed in models of turbulent ad- 
vection, 

We know however, that mean field approximation is 
invalid in the limit of large times if dimension is two or 



less because of strong fluctuations of local particle density 
. Using results of the previous sections one can interpret 
Kang-Redner anomaly as the anomaly in the constant 
flux condition: 



dJ_ 



J, 



(50) 



where ckr is Kang-Redner's exponent. 

Solving Eq. (|50|), we find that P{m, t) ~ /x*^, where e = 
Skoim + Therefore, Kang-Redner anomaly can be 

also interpreted as a correction to Kolmogorov scaling of 
the average mass distribution due to strong fluctuations. 



VII. CONCLUSION. 

In the present paper we have shown that the problem 
of cluster-cluster aggregation in d < 2 can be effectively 
analysed using renormalization group method. We have 
demonstrated that the dependence of average mass dis- 
tribution on mass is determined by the anomalous di- 
mension of the stochastic field P (local mass distribu- 
tion). This anomaly is due to the relevance of 'bound- 
ary' {t = 0) fluctuations for the large times asymptotics 
of P(m,t). In that respect the phenomenon of Kang- 
Redner anomaly resembles the phenomenon of boundary 
phase transition in equilibrium statistical mechanics, see 
Q for a review. Formally, the anomalous dimension of 
the local mass distribution is a consequence of the non- 
triviality of the 7-function of the effective field theory 
Eq. (p^. The fact that 7(5) ^ is ultimately respon- 
sible for the breakdown of the Smoluchowski theory (or 
equivalently, the renormalized mean field theory) applied 
to the model at hand. At present, renormalization group 
analysis seems to be the only theoretical method of study- 
ing the problem of cluster-cluster aggregation in d > 1. 
This is not quite satisfactory, as the analysis is essentially 
perturbative in nature. However, our theoretical predic- 
tions concerning the behavior of the average mass dis- 
tributions at small masses in two dimensions have been 
unambiguously confirmed numerically. 

have shown 



In section VI 



we have stiown that there is a rela- 
tion between cluster-cluster aggregation and the theory 
of weak turbulence. In particular, we demonstrated that 
Kang-Redner anomaly can be interpreted as a correc- 
tion to the Kolmogorov spectrum of particles in the mass 
space. However, we must stress that the analogy between 
our model and the phenomenon of turbulence must be 
taken with a pinch of salt: while the cascade of the con- 
served quantity in our model happens along the mass axis 
only, conserved quantities (such as energy or enstrophy) 
in more traditional turbulent systems flow through the 
scales of the physical space. 

The method of dynamical renormalization group de- 
veloped in the context of the model at hand can be ap- 
plied to other non-equilibrium particle systems as well. 
In particular the problem of cluster-cluster aggregation 
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with annihilation of particles, can be solved using tech- 
niques similar to those outlined in this paper pl| . This 
latter problem is related to the computation of the do- 
main wall persistence exponent for the Id q-state Potts 
model g. 
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